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Charmonia spectral functions at finite temperature are studied using QCD sum rules in combina- 
tion with the maximum entropy method. This approach enables us to directly obtain the spectral 
function from the sum rules, without having to introduce any specific assumption about its func- 
tional form. As a result, it is found that while J/tp and rjc manifest themselves as significant peaks 
in the spectral function below the deconfinement temperature Tc, they quickly dissolve into the 
continuum and almost completely disappear at temperatures between 1.0 Tc and 1.1 Tc. 
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Since QCD was established to be the theory of strong 
interactions, charmonium has often been used as a suit- 
able probe of its dynamics, owing to the fact that in this 
system both pcrturbative and nonperturbative aspects of 
QCD play equahy important roles [l|. The behavior of 
charmonia in a hot or dense medium has also attracted 
much interest, as it was suggested some time ago, that in 
the color-deconfined medium with a temperature above 
Tc charmonia will dissolve due to the color Debye screen- 
ing, and thus serve as a signal for the formation of quark- 
gluon plasma Q ■ 

Testing these early suggestions from first principles of 
QCD has become feasible only recently, as new devel- 
opments in lattice QCD have made it possible to access 
the charmonium spectral functions with the help of the 
maximum entropy method (MEM) [3^^. These studies 
found that the lowest charmonium states {J/ip and rjc) 
survive up to temperatures as high as ^ 1.5 Tc or higher. 

Besides lattice QCD, the method of QCD sum rules 
provides another tool for investigating the proper- 
ties of hadrons at finite temperature [1, Q. Using this 
approach various charmonium channels were studied re- 
cently [13, 11|, and evidence for a considerable change 
of the spectral functions just above Tc was found. To 
specify the nature of this change is the major goal of this 
study. For this task we employ MEM, which is applicable 
to QCD sum rules [12] and has the advantage that one 
does not have to introduce any strong assumption about 
the functional form of the spectral function, such as the 
"pole + continuum" ansatz, which is often used in QCD 
sum rule studies. 

Let us first recapitulate what sort of information QCD 
sum rules can provide on the charmonium spectral func- 
tion at finite temperature [^, . One considers the time- 
ordered correlator at finite temperature 

U\q) = I J d^xe^'^^T[j\x)j\0)])T, (1) 



where j'^{x) stands for 07^0(2;) and C75c(a;) in the vec- 
tor (V) and pseudoscalar (PS) channel, respectively. 
The expectation value {0)t is defined as {0)t = 
TT{e-"/'^0)/TT{e-"/'^). Throughout this work, we will 
set the spatial momentum of the charmonium system rel- 
ative to the thermal medium to be 0; thus, = {^,0)- 
In this circumstance, there is only one independent com- 
ponent in the correlator of the vector channel. In what 
follows, we will use the dimensionless functions 11^ (q^) = 
n;;-^(g)/(-3g2) and Il^^{q^) = Il^^{q)/q^ for the analy- 
sis. 

Going to the deep Euclidean region q'^ = —Q^ -C 0, 
one can calculate the correlation functions using the op- 
erator product expansion (OPE), giving an expansion in 
local operators 0„ with increasing mass dimension n: 
IlJ(q2) ^ Ci{q^){On)T- As was first discussed in [3, 
as long as the temperature T lies below the separation 
scale of the OPE, which is of the order of 1 GcV, all 
the temperature effects can be included into the expecta- 
tion values of the local operators {On)T, while the Wilson 
coefficients C:l{q^) are independent of T. Furthermore, 
to improve the convergence of the OPE and suppressing 
the infiuence of high energy states onto the sum rule, we 
apply the Borel transform to the correlator, leading to 
the final result of the OPE for = 4ml /M^, ^ being 
the Borel mass: 



+ c\u)cl,c{T) + d\v)<j,d{T)]. 



(2) 



The first two terms in Eq. ^ are the leading order pcrtur- 
bative term and its first order as correction. The third 
and fourth terms contain the scalar and twist-2 gluon 

0h(r) = 37^7^ Go 



condensates of mass dimension 4 
and (pciT) - 



9(4m^)^ ' 

j^G2, where Go = {^Gl^G^nr and 

G2 is defined as {^G''^"'G%'')T = (u'^u'' - iff'''')G2, u^' 
being the four velocity of the medium. For the detailed 
expressions of the Wilson coefficients of these terms, see 
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To evaluate the possible influence of higher order 
contributions, we include one more term, which is pro- 
portional to the scalar gluon condensate of dimension 6, 
UT) = (iT^CgV^'^G^'-G^^Gf )t. The Wilson coeffi- 
cient of this term can be found in 



13|. 



The correlator can also be expressed by a dispersion 
relation, in terms of the spectral function of the 

channel specified by the operator j'^{x). After the Borel 
transform one obtains 



(3) 



Equating Eqs.® and ^ then gives the final form of the 
sum rules. In the vector channel, an additional constant 
term contributes to Eq.([5]), which originates from a pole 
at w = in p^{i-u) As this so-called scattering term 
considerably complicates the analysis, we eliminate it by 
taking the derivative of Eqs.® and ([3]) with respect to 
and analyze only the resulting derivative sum rule in this 
channel. For a discussion on the validity of this procedure 
in the heavy quark sum rules, see [l^ . 

The usual strategy of analyzing QCD sum rules is 
to make some reasonable assumptions on the functional 
form of the spectral function, and then extract informa- 
tion on the lowest lying peak from Eqs.(l2]) and ((3)). This 
method, however, has several shortcomings. First of all, 
the widely used "pole -I- continuum" ansatz, which cer- 
tainly works well at T = 0, may not be appropriate at 
temperatures above Tc, where the lowest lying state is 
expected to be modified and eventually melt into the c-c 
continuum, which could become the dominant contribu- 
tion. Furthermore, it is not always possible to unam- 
biguously fit a specific ansatz to the OPE results, be- 
cause of the occurrence of equally valid solutions. Such 
a situation arose in lol 11| 



where it was not possible 
to determine a unique solution for the used parametriza- 
tion of the spectral function. To handle these problems, 
we propose to use MEM, which allows us to extract the 
spectral function from Eqs.® and @ without prejudice 
on its functional form. Moreover, it can be proven that 
this method provides a unique solution for the spectral 
function 
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Let us now briefly summarize the basic ideas of MEM, 
which helps us to carry out the task of inverting the in- 
tegral of Eq.(j31). This is, however, an ill-posed problem 
as we have only access to A4'^{i') in the region of ly where 
the OPE shows sufficient convergence and, furthermore, 
have only information on Ai'^{i') with limited precision 
due to the uncertainties of the values of nic, ots and the 
various condensates. Nonetheless, using Bayesian prob- 
ability theory, MEM makes it possible to at least obtain 
the most probable form of the spectral function piuj). 
To do this, one defines a probability P[p|A4I] for p to 
have a specific form given Ai and additional information 
I on p such as positivity and asymptotic values. Us- 
ing Bayes' theorem, this probability can be rewritten as 



P[p\MX] oc P[M\pI]P[p\I] = e-^+"'S', where e'^ is the 
likelihood function, used in standard methods, e"'^ 
stands for the prior probability, given by the parameter 
a and the Shannon- Jaynes entropy, 



S 



du 



pH-™H-pHiog(^) 



(4) 



m{uj) is called the default model and can be used to in- 
corporate prior knowledge on the asymptotic values of 
p{uj) into the analysis. Determining now the most prob- 
able p{ll!) corresponds to the numerical problem of finding 
the maximum of the functional —L + aS, for which we 
use the Bryan algorithm [l^. As a last step, the spectral 
function Pa{<^) maximizing —L-\-aS for a fixed value of 

a is averaged using j c?apQ,(aj)P[Q!|A^I], which gives the 

final output of the MEM analysis. For more details see 
15l |. while specific issues concerning the application of 
this method to QCD sum rules are discussed in jl2|. 

Next, we describe how the temperature dependencies 
of the gluonic condensates are determined. For the scalar 
and twist- 2 gluon condensates with mass dimension 4, 
we follow the approach proposed in 0,1 where, in 
the quenched approximation, the energy momentum ten- 
sor, expressed using gluonic operators, was matched with 
the corresponding quantity written down in form of the 
energy density e and the pressure p, leading to Go = 



-A[e(T)™3p(T)] andG2 



s(T) 



'e{T)+p{T)\ for 



Gvac 
' 

the scalar and twist- 2 gluon condensates. The functions 
e(T), p{T) and a^fr) were then extracted from quenched 
lattice QCD data [11 El. We will in this study use the 
same numerical values for the T dependent part of Gg 
and G2 as in As is shown there, both Go and G2 

exhibit a sudden decrease in the vicinity of Tc- 

It was suggested in previous studies that the OPE 
could break down at temperatures above Tc as higher di- 
mensional operators may become non- negligible To 
investigate this possibility, we include the scalar gluonic 
condensate with dimension 6, {g^ /"'""C^'^Gl^Gl^), about 
which much less is known. To our knowledge, at T = 0, 
there exists only an estimate based on the dilute instan- 
ton gas model, giving {g"^ f^^C^^Gl^Gf) = f^{^G^), 
where Pc is a representative value for the instanton ra- 
dius, for which we use the established value of 0.33 fm 
[l9| . Assuming that the relation above also holds at finite 
temperature, and taking into account the reduction of pc 
above Tc [20[, we, however, found that the dimension 6 
term does not influence the behavior of the spectral func- 
tion much in the temperature region investigated in this 
Letter. Therefore, we conclude that even though the re- 
lation obtained from the dilute instanton gas model can 
only be considered to be a crude estimate, as long as it 
gives the correct order of magnitude, the contribution of 
the dimension 6 condensate is small and does not alter 
the results obtained in this study. 
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FIG. 1: On the left side, the spectral functions in the pseudoscalar (upper plot) and vector (lower plot) channel at T = 0, 
with the errors of their averaged values are shown. The horizontal extent and position of the error bars indicates the region 
over which the average is taken. The dashed lines show the default model used in the MEM analysis. On the right side, the 
corresponding spectral functions are given at temperatures around Tc. 



Let us now turn to the MEM analysis of Eqs.(I2|) and 
First, wc investigate the spectral function at T = 
in the vector and pseudoscalar channel. To determine 
the upper boundary of the region of v to be analyzed, we 
employ the criterion that the dimension 6 term should be 
smaller than 20% of the whole OPE expression of Eq. ^ , 
which gives f^ax = 8-03 in the vector and fJnlx = 7-29 
in the pseudoscalar channel. We keep these values fixed 
when going to finite temperature. In fact, in the tempera- 
ture region around Tc, the relative contribution of the di- 
mension 6 term at t'^ax^ is even smaller, namely, around 
10% or less. The lower boundary of v is chosen to be 

V PS 

'^min = 0.78, corresponding to a Borel mass of M — 3.0 
GeV. We have checked that the obtained spectral func- 
tions do not depend much on this choice. For the value of 
the charm quark mass m^, we use a recent estimate giv- 
ing mc(mc) = 1.277 ± 0.026 GeV [lH, for we emplw 
the newest world average as{Mz) = 0.1184±0.0007 |U, 
while for the vacuum gluon condensate G™"^ the standard 
value G^^'= = 0.012 ± 0.0036 GeV"^ 0,|23l is applied. For 
the default model m(w), we use a constant matched to 
the perturbative value of the spectral function at high 
energy, as was done in similar studies using lattice QCD 

i- 

The resulting spectral functions are given on the left 
side of Fig. [TJ We observe in both channels a clear ground 
state peak, corresponding to rjc and J/V'- The spectral 
functions also exhibit a second peak, which is, however, 
not statistically significant. These second peaks most 
likely reflect the existence of several excited states, which 
the MEM analysis is not able to resolve, quite similar to 



the situation encountered in lattice studies. Furthermore, 
it is seen that the spectral function of the vector channel 
approaches the default model in the region close to w = 0, 
which; however, should not be confused with a contribu- 
tion of the scattering term. This behavior is an artifact 
caused by our usage of the derivative sum rule in this 
channel and should thus not be considered to be a phys- 
ical effect. Numerically, the peak representing 77c lies at 
3.02 GeV, while the one standing for J/ip is found at 3.06 
GeV. Thus, we see that the ground state in both channels 
reproduces the experimental value with a precision of the 
order of 50 MeV. In the vector channel, the residue can 
be related to the electronic width of the corresponding 
resonance. We can obtain this residue from Fig. [T]simply 
by integrating the spectral function in the region of the 
peak, which gives 0.162 GeV^. which is in good agree- 
ment with the experimental value of 0.173 GeV^. On the 
other hand, we observe that the hyperfine splitting be- 
tween rjc and J /ip is underestimated. All these findings 
are in qualitative agreement with the results obtained in 
the conventional analysis of the charmonium sum rules 

Next, we increase the temperature according to Eq.([2|). 
The resulting spectral functions are shown on the right 
side of Fig. [T]at temperatures between 0.9 Tc and 1.2 Tc. 
It is seen in the figure that the behavior of the spectral 
functions changes abruptly in the vicinity of T^. First, 
both ground state peaks experience a shift to lower en- 
ergies of the order of 50 MeV, before dissolving quickly 
into the continuum above the critical temperature. In- 
vestigating the spectral functions in more detail, one ob- 
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serves that rjc disappears already at T = 1.0 Tc, while 
J/ip survives a bit longer, but also appears to be melted 
to a large degree at T = 1.1 T^. This sudden quali- 
tative change of the spectral function mainly originates 
from the changes of the third and fourth terms in Eq. 1^ , 
which can be traced back to the rapid adjustment of the 
thermodynamic quantities e{T) and p{T) around Tc- It 
is reassuring to note that our results are consistent with 
the findings of 13, III in the sense that both observe a 
negative energy shift of the peaks around Tc- In these 
earlier works, it was, however, not possible to discuss the 
possible melting of the peaks because a relativistic Breit- 
Wigncr form for the spectral function was assumed at all 
investigated temperatures. 

For obtaining firm conclusions, one has to test the re- 
liability of the MEM procedure at finite temperature, 
where systematic effects decrease the reproducibility and 
resolution of the spectral function obtained from MEM. 
In lattice studies, this reduced reliability is primarily 
caused by the reduction of the data points in the imag- 
inary time correlator, due to periodicity and the reduc- 
tion of the maximal time extent. In the case of QCD sum 
rules, this problem does not exist, as Eq.([2]) is given as 
a continuous function at any temperature and therefore 
the same number of data points can be used. Neverthe- 
less, the reliability of the MEM technique is still reduced 
at finite temperature due to the uncertainties of the ther- 
modynamic functions e{T) and p(T), whose contribution 
grows with temperature and therefore increases the er- 
ror of Eq. ([2]) . In order to confirm that the change of the 
spectral function in Fig. [T] is not an artifact, we rean- 
alyze Eq.© at T = 0, but use the errors of T 7^ in 
the analysis. The results are then compared to the ones 
given in Fig. [1] to investigate the net temperature effect 
on the spectral function. We find from this analysis that 
while the height of the peaks of the spectral functions at 
T = is indeed reduced because of the increased error, 
this effect is much smaller than the actual reduction of 
the peaks around Tc, seen on the right side of Fig. [TJ We 
therefore conclude that the disappearance of the peaks 
observed in Fig. [1] is a physical effect and is not induced 
by an artifact of the MEM analysis. 

In summary, we have extracted the spectral functions 
of the pseudoscalar and vector channel at both zero and 
finite temperature using a combined analysis of QCD sum 
rules and MEM. At T = 0, the MEM technique is able to 
clearly resolve the lowest energy peaks, corresponding to 
the rjc and J/^p resonances. The positions of both peaks 
agree with the experimental values with a precision of 
about 50 MeV. At finite temperature, we find that ijc and 
J/ip melt quickly after the temperature is raised above 
the deconfinement temperature Tc, caused by the sudden 
change of the dimension-4, scalar and twist- 2 gluon con- 
densates in this temperature region. We have checked 
that this effect is not an artifact of the systematics of 



the MEM analysis. These results quantitatively disagree 
with the earlier findings of lattice studies which suggest 
that both ric and J/^p can survive at temperatures of up 
to 1.5 Tc or higher. It, however, has to be mentioned 
that our results are in fact consistent with the latest lat- 
tice results [i^l, finding the peaks of rjc and J/^j to be 
largely distorted between 0.73 Tc and 1.46 Tc- It remains 
to be seen whether or not the two methods will converge 
to compatible conclusions when more accurate analyses 
will become available in the future. 
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